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Abstract 

In this article, we present the analytical solution of the radial Schrodinger equation for the 
Hulthen potential within the framework of the asymptotic iteration method by using an approx- 
imation to the centrifugal potential for any I states. We obtain the energy eigenvalues and the 
corresponding eigenfunctions for different screening parameters. The wave functions are physical 
and energy eigenvalues are in good agreement with the results obtained by other methods for dif- 
ferent 5 values. In order to demonstrate this, the results of the asymptotic iteration method are 
compared with the results of the supersymmetry, the numerical integration, the variational and 
the shifted 1/N expansion methods. 
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I. INTRODUCTION 



Over the last decades, the energy eigenvalues and corresponding eigenf unctions between 
interaction systems have raised a great deal of interest in relativistic quantum mechanics 
as well as in non-relativistic quantum mechanics. The exact solution of the wave equa- 
tions (relativistic or non-relativistic) are very important since the wave function contains all 
the necessary information regarding the quantum system under consideration. Analytical 
methods such as the supersymmetry (SUSY) [lj and the Nikiforov-Uvarov method (NU) [J| 
have been used to solve the wave equations exactly or quasi-exactly for non zero angular 
momentum quantum number (I ^ 0) by means of a given potential. The radial Schrodinger 

nnn 

equation for the Hulthen potential is solved exactly by using several techniques ja 121 ISI] for 

1 = 0. For the case / 7^ 0, the effective Hulthen potential can not be solved exactly, but a 

number of methods have been used to find the boiind state energy eig envalues numerically 

perturbation lOj, shifted 1/N expan- 
, , ' 

sion 



and quasi-analytically such as the variational 

NU and SUSY methods. The Hulthen potential [2] is one of the important 
short-range potentials in physics and it potential has been applied to a number of areas such 
as nuclear and particle physics, atomic physics, condensed matter and chemical physics (see 
and the references therein). Therefore, it would be interesting and important to solve 
the non-relativistic radial Schrodinger equation for this potential for Z 7^ 0, since it has been 
extensively used to describe the bound and the continuum states of the interaction systems. 
Recently, an alternative method, called as the asymptotic iteration method (AIM), has been 
developed by Ciftgi et al. ^, "\ for solving second-order homogeneous linear differential 
equations and it has been applied to solve the non-relativistic radial Schrodinger equation 
or the relativistic Dirac equation. 

In this paper, we aim to solve the Hulthen potential to obtain the energy eigenvalues 
and corresponding eigenf unctions for any I states. In the next section, AIM is introduced. 
Then, in section the Schrodinger equation is solved by using AIM for the Hulthen 
potential for any I states and our AIM results are given in comparison with the results of 
the numerical integration [9], the variational the shifted 1/N expansion and the 
SUSY [3] methods. Finally, section |IV| is devoted to the summary and conclusion. 
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II. OVERVIEW OF THE ASYMPTOTIC ITERATION METHOD 



A. Energy Eigenvalues 

AIM is briefly outlined here and the details can be found in references P, 1^, |5|. AIM is 
proposed to solve the second-order differential equations of the form 

y" = X (x)y' + s (x)y (1) 

where \o(x) ^ and the prime denotes the derivative with respect to x. The variables, 
s (x) and_^Ao(x), are sufficiently different iable. The differential equation (JTJ) has a general 
solution 



^— J a^x^dx-^j C2 + C1 J exp(^J [Xq(x2) + 2a(x2)]dx2*j dx\ 



y(x) = exp ^— J a(xi)dxi 

for sufficiently large k, k > 0, if 

Sk{x) s fe _i(x) 



(2) 



X k (x) Xk-i(x) 



a(x), k = 1,2,3,... (3) 



where 



A fc (x) = X'^ix) + s fe _i(a;) + A (a;)A fc _i(a;) 

s k(x) = 4-i C^) + s (x)X k - 1 (x), fc = 1,2,3,... (4) 

Note that one can also start the recurrence relations from k = with the initial conditions 
A_i = 1 and s_! = j^J. For a given potential such as the Hulthen, the radial Schrodinger 
equation is converted to the form of equation (0). Then, sq(x) and Xq(x) are determined and 
Sfc(x) and Afc(x) parameters are calculated by the recurrence relations given by equation 

The termination condition of the method in equation (jSJ) can be arranged as 

Afc(x) = A fc (x)s fe _i(a;) - A fc _i(x)s fc (x) = k = 1,2,3,... (5) 

The energy eigenvalues are obtained from the roots of the equation (j5J) if the problem is 
exactly solvable. If not, for a specific n principal quantum number, we choose a suitable Xq 
point, determined generally as the maximum value of the asymptotic wave function or the 
minimum value of the potential 0,13, 19, an d the approximate energy eigenvalues are 



obtained from the roots of this equation for sufficiently great values of k with iteration. 



B. Energy Eigenfunctions 



In this study, we seek the exact solution of the radial Schrodinger equation for which the 
relevant second order homogenous linear differential equation takes the following general 



form 4|, 



,. , , ax"* 1 (t + m , w' k (N)x N 



If this equation is compared to equation ((TJ, it entails the following expressions 

, / X nf axN+1 ( t + 1 )\ , N W{{N)X N 

a and 6 are constants and w\{N) can be determined from the condition (jSJ) for k = 0, 1, 2, 3, ... 



and iV=-l,0,l,2,3,... as follows: 



4(-i) 


= fc (2a + 26t + (jfc + 1)6) 


4(o) 


= 2k (2a + 2bt + (2k + 1)6) 




= 3k (2a + 2bt + (3k + 1)6) 


4(2) 


= 4k (2a + 2bt + (4k + 1)6) 


4(3) 


= 5k (2a + 2bt + (5k + 1)6) 


... etc 





(8) 
(9) 
(10) 

(11) 
(12) 



Hence, these formule are easily generalized as, 



m < m = H N + 2fk(k + ^0^ ) (13) 
The exact eigenfunctions can be derived from the following generator: 

y n (x) = C 2 exp f ^^ dx ^j ( M ) 

where k >n,n represents the radial quantum number and k shows the iteration number. For 
exactly solvable potentials, the radial quantum number n is equal to the iteration number 
k and the eigenfunctions are obtained directly from equation (J14J) . For nontrivial potentials 
that have no exact solutions, k is always greater than n in these numerical solutions and the 
approximate energy eigenvalues are obtained from the roots of equation (jSJ) for sufficiently 
great values of k with iteration. It should be pointed out that a(x) given by equation (j3J) is 



equal to zero for the ground state. Therefore, using equation (jSJ) with (JJJ) in equation (JHJ), 
the eigenf unctions are obtained as follows, 



y (x) = G 



2 



Vi{x) = -C 2 (N + 2)a[l- ^ > x 



b (p + 1 ) „N+2 



V °" cr(o-+l) J 



( v a (<7 + 1) (a + 2 

2/3 W = -C 2 



(JV + 2)~ 3 

3fe(p + 3) „ +a 36 2 (p + 3) (p + 4) 2(JV+2) _ fe 3 (p + 3) (p + 4) (p + 5) 3(JV+2 A 
' a X + <r(a + l) X a(a + l)(a + 2) X J 

. . . etc 

Finally, the following general formula for the exact solutions y n (x) is obtained as, 

y n (x) = (-l) n C 2 (N + 2)" {a) n 2 F 1 (-n, p + n;a; bx N+2 ) (15) 

It is important to note that square integrable one in L 2 is this total wave function which 
is the asymptotic form of the wave function times y n (x) given by equation (jl4j) . Here, 

(°")n— re?) i ' a = 2t jv+2~ 3 an< ^ <° = ^JW+^F^- ^ e (° r ) n an< ^ ^he 2 -^ 1 are known as the 
Pochhammer symbol and the Gauss hypergeometric function, respectively. 



III. CALCULATION OF THE ENERGY EIGENVALUES AND EIGENFUNC- 
TIONS 

The motion of a particle with the mass M in the spherically symmetric potential is 
described in the spherical coordinates by the following Schrodinger equation: 



-h 2 ( d 2 2 8 1 



2M V dr 2 r dr r 2 



1 d 
sin 9 89 



sin9-^- ] + 



1 d 2 



V(r) ) ^ nlm (r,9,(f)) = E^ nlm (r } 9 } 

(16) 



Defining ^ nim (r, 6>, 0) = R n i(r)Y lm (9, 0), 
tion: 

d 2 2 d \ 2M 
dr 2 r dr J nl h 2 



89 J sm 2 9d(f) 2 
we obtain the radial part of the Schrodinger equa- 



E-V(r) - 



l(l + i)h 2 

2Mr 2 



Rm(r) = 



(17) 



5 



It is sometimes convenient to define Rni(r) and the effective potential as follows 

«-(D = ^, vv,, = vM + !2±a* (18 ) 

Since 

d 2 2 d \ Uniir) 1 d 2 



dr 2 ^ r dr J r r dr' 2 ^ ^ ^ 



The radial Schrodinger equation |8( given by equation (|T7|) follows that 

d 2 u nl (r) , 2M 



dr 2 ft 2 
The Hulthcu potential (1 is given by 



+ — [E-V eff ]u nl (r) = (20) 



e -8r 



V H (r) = -Ze^^—^ (21) 

where Z and 5 are respectively the atomic number and the screening parameter, determin- 
ing the range for the Hulthen potential. The Hulthen potential behaves like the Coulomb 
potential near the origin (r — > 0), but in the asymptotic region (r ^> 1), the Hulthen po- 
tential decreases exponentially, so its capacity for bound states is smaller than the Coulomb 
potential. However, for small values of the screening parameter 5, the Hulthen potential 
becomes the Coulomb potential given by Vc = The effective Hulthen potential is 



K,,(r) = V„(r) + V, = -Ze'Sj^ + '-^±±f (22) 

where Vi = '^ 2 fi is known as the centrifugal term. This effective potential can not be 
solved analytically for I ^ because of the centrifugal term. Therefore, we must use an 
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13 



14. 



Q0,[i2|. 



In 



approximation for the centrifugal term similar to other authors 
this approximation, ^ = 5 2 ^ 1 _ e e _L)2 is used for the centrifugal term. As shown in figure [U 
this approximation is only valid for small 5r and it breaks down in the high-screening region. 
For small 8r, V e ff(r) is very well approximated to V e ff(r) and the Schrodinger equation for 
this approximate potential is solvable analytically. So, the effective potential becomes 



V eff (r) = -Ze 2 5—— + ^^ L > ,bU (23) 



e- Sr 1(1 + l)h 2 5 2 e~ 6r 
1 - e~ Sr + 2M (1 -e" 
Instead of solving the radial Schrodinger equation for the effective Hulthen potential 
V e ff(r) given by equation (|22|). we now solve the radial Schrodinger equation for the new 
effective potential K//( r ) given by equation (J2BJ). Inserting this new effective potential into 



equation (j20j) and using the following ansatzs in order to make the differential equation more 
compact 



2ME 



2 



2MZe 2 



h 2 5 2 ' y h 2 5 ' 
The radial Schrodinger equation takes the following form: 

d 2 u n i(x) 



5r = x 



dx 2 



1(1 + 1) 



u ni (x) = 



1 - e~ x ) " v ' ' (1 - e-*) 2 _ 
If we rewrite equation (j2*Hj) by using a new variable of the form z = e~ x , we obtain 

2 (3 2 1(1 + 1) 



(24) 



(25) 



d 2 u n i(z) lduni(z^ 



dz 2 



z dz 



e 

— o + 



z(l-z) z(l-2) 



u n i(z) = 



(26) 



In order to solve this equation with AIM, we should transform this equation to the form of 
equation (JTJ). Therefore, the reasonable physical wave function we propose is as follows 



u nl (z) = z £ (l-z) l+l f nl (z) 



(27) 



If we insert this wave function into equation 1)26)1 . we have the second-order homogeneous 
linear differential equations as in the following form 



d 2 fnl(z) 

dz 2 



(2e + 21 + 3)z - (2e + 1) 



z(l-z) 



dfnl(z) 

dz 



+ 



(2s + 1 + 2)1 + 2e-(3 2 + l 



z(l-z) 



fm(z) (28) 



which is now amenable to an AIM solution. By comparing this equation with equation ( ) T ]) . 
we can write the A (z) and s (z) values and by means of equation (JJJ), we may calculate 
Afc(z) and Sk(z). This gives: 

. . , ((2e + 2l + 'i)z-(2e+\) 
X 1 (z) 

Z" {-L + Z) 

12^2 



z(l - z) 
(2e + l + 2)1 + 2e - f3 2 + 1 

z(l ~ *) 

2 + 6 e - 7 z - 2 Iz - /3 2 z + 12 z 2 l - 18 e z - 6 e zl 
z 2 (-l + z) 2 



Sl (z) 
. . . etc 



12 e z 2 + 11 z 1 + 4 e 1 + Vz + (3 2 z 2 + 4 e V - 8 e 2 z + 6 e zH + 3 Vz 

z 2 (-\ + zf 

(21 + 2e - (3 2 + 2e I + I 2 + I) (-2 + 5 z + 2e z + 2lz - 2e) 

z 2 (-l +zf 



(29) 
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Combining these results with the quantization condition given by equation (jSJ) yields 

(3 2 - 1 - 21 - I 2 



s Ai - s 1 X = 

siA 2 - s 2 Ai = 

S2A3 — S3A2 = 
. . . etc 



£0 



£1 



£2 



2(1 + 1) ' 

2(1 + 2) 
f3 2 -9-Ql-l 2 
2(1 + 3) : 



for k = 1 
for k = 2 
for k = 3 



(30) 



When the above expressions are generalized, the eigenvalues turn out as 

'(3 2 -(n + l + l) 2 * 



£nl 



2(n + l + l) 

Using equation ()24|1. we obtain the energy eigenvalues E n/ 



n,l = 0,1,2,3, 



(31) 



E„ 



h 2 



nl 



2M 



MZe 2 



(n + 1 + 1)5 



h 2 (n + l + l) 



(32) 



In the atomic units (h — M — e — 1) and for Z = 1, equation (p?2"|) turns out to be 



E 



nl 



1 



(n + l + l)8 



(33) 



(n + l + 1) 2 

In order to test the accuracy of equation (J3~3"j) . we calculate the energy eigenvalues for 
Z — 1, any n and I quantum numbers and several values of the screening parameter. AIM 
results are compared with the results of the numerical integration |QJ, the variational 
the shifted 1/N expansion 11] and the SUSY jl^ methods in Table |J and Table ITT1 As it 
can be seen from the results presented in these tables, the AIM results are in good agrement 
with the results of the other methods for the small 5 values. For large 5 values, there are 
differences between our results and the results of others. This difference is due to the V e ff(r) 
potential, which we have used to approximate the V e f/(r) potential. As it is seen from figure 
[TJ for large 5r values, the discrepancy becomes apparent between our V e ff(r) and the true 
V e ff(r) potentials. This gives rise to the differences for the eigenvalues presented in Tables 
HI and ITT1 at large 6 values. 

Now, as indicated in Section |Hj we can determine the corresponding wave functions by 
using equation ([To]) . When we compare equation (0) and equation (|2"H|) . we find N = —1, 
b = 1, a = I + 1, and t = Therefore, we find p = 2(e + 1+1) and a = 2e + 1. So, 

we can easily find the solution for f n i(z), for the energy eigenvalue equation (JH2~|) by using 
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equation (|T5|). 

f nl {z) = (-l)" r( ^" + ^+ 1) 2 F 1 (-n, 2e„ + 2/ + 2 + n; 2e n + 1; z) (34) 
Thus, we can write the total radial wave function as below, 

= Nz £ - (1 - z)' +1 2 Fi(-n, 2(e n + I + 1) + n; 2e n + 1; z) (35) 
where A/" is the normalization constant. 

IV. CONCLUSION 

We have shown an alternative method to obtain the energy eigenvalues and corresponding 
eigenfunctions of the Hulthen potential within the framework of the asymptotic iteration 
method for any / states. We have calculated the energy eigenvalues for the Hulthen potential 
with Z = 1 and several values of the screening parameter. The wave functions are physical 
and energy eigenvalues are in good agreement with the results obtained by other methods. In 
order to demonstrate this, AIM results have been compared with the results of the numerical 
integration |9j , the variational , the shifted 1/N expansion Q and the SUSY methods 
in Tables U and |H] For small 5 values, AIM results are in good agreement with the results 
of the other methods, but in the high screening region, the agreement is poor. The reason 
is simply that when the Sr increases in the high screening region, the agreement between 
V e ff(r) and V e ff(r) potentials decreases as shown in Figure HJ This problem could be solved 
by making a better approximation of the centrifugal term. 

It should be pointed out that the asymptotic iteration method gives the eigenvalues 
directly by transforming the radial Schrodinger equation into a form of y" =Ao(?")?/ + so(r)y. 
The wave functions are easily constructed by iterating the values of sq(t) and Xo(r). The 
asymptotic iteration method results in exact analytical solutions if there is and provides 
the closed-forms for the energy eigenvalues as well as the corresponding eigenfunctions. 
Where there is no such a solution, the energy eigenvalues are obtained by using an iterative 
approach As it is presented, AIM puts no constraint on the potential parameter 



values involved and it is easy to implement. The results are sufficiently accurate for practical 
purposes. 
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State 


<5 


AIM 


SUSY [13] 


Numerical Int. [9] 


Variational [9] 


Shifted i [111 

iv ^« 


2p 


0.025 


0.1128125 


0.1127605 


0.1127605 


0.1127605 






0.050 


0.1012500 


0.1010425 


0.1010425 


0.1010425 


0.1010424 




0.075 


0.0903125 


0.0898478 


0.0898478 


0.0898478 






0.100 


0.0800000 


0.0791794 


0.0791794 


0.0791794 


0.0791794 




0.150 


0.0612500 


0.0594415 


0.0594415 


0.0594415 






0.200 


0.0450000 


0.0418854 


0.0418860 


0.0418860 


0.0418857 




0.250 


0.0312500 


0.0266060 


0.0266111 


0.0266108 






0.300 


0.0200000 


0.0137596 


0.0137900 


0.0137878 






0.350 


0.0112500 


0.0036146 


0.0037931 


0.0037734 




3p 


0.025 


0.0437590 


0.0437068 


0.0437069 


0.0437069 






0.050 


0.0333681 


0.0331632 


0.0331645 


0.0331645 


0.03316518 




0.075 


0.0243837 


0.0239331 


0.0239397 


0.0239397 






0.100 


0.0168056 


0.0160326 


0.0160537 


0.0160537 


0.01606772 




0.150 


0.0058681 


0.0043599 


0.0044663 


0.0044660 




3d 


0.025 


0.0437587 


0.0436030 


0.0436030 


0.0436030 






0.050 


0.0333681 


0.0327532 


0.0327532 


0.0327532 


0.0327532 




0.075 


0.0243837 


0.0230306 


0.0230307 


0.0230307 






0.100 


0.0168055 


0.0144832 


0.0144842 


0.0144842 


0.0144842 




0.150 


0.0058681 


0.0132820 


0.0013966 


0.0013894 





TABLE I: Energy eigenvalues of the Hulthen potential as a function of the screening parameter 
for 2p, 3p and 3d states in atomic units (h = m = e = 1) and for Z = 1. 
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State 


s 


AIM 


SUSY [131 


Numerical Int. [9] 


Vaxia/tioiicil [9] 


Shifted 4? [Ill 


4r> 


025 


0200000 


01 99480 


01 QQ48Q 


D 01QQ48Q 








01 1 2500 


01 1 04S0 


01 1 0582 


01 1 0582 


01 1 0725 




0.075 


0050000 


0045385 


0046219 


0046219 






0.100 


0012500 


0004434 


0007550 


0007532 




4d 


0.025 


0200000 


0198460 


0198462 


0198462 






0.050 


0112500 


0106609 


0106674 


0106674 


0106690 




0.075 


0.0050000 


0.0037916 


0.0038345 


0.0038344 




4f 


0.025 


0.0200000 


0.0196911 


0.0196911 


0.0196911 






050 


01 1 2500 


01 0061 8 


01 00620 


01 00620 


01 00620 




0.075 


0050000 


0025468 

\J •\J\J £iKJ^.\JkJj 


0025563 


0025557 




5n 
°i' 


0.025 


0094531 

\J •\J\J tJ^tJ\J ±. 


0094011 

\J •\J\J U^.\J ±- X 


0094036 




0094087 




0.050 


0028125 


0026056 


0026490 






5d 


0.025 


0094531 


0092977 


0093037 




0093050 




0.050 


0028125 


0022044 


0023131 






5f 


0.025 


0.0094531 


0.0091507 


0.0091521 




0.0091523 




0.050 


0.0028125 


0.0017421 


0.0017835 








0.025 


0094531 

\J •\J\J U^.tJKJ -L 


0089465 


0089465 




0089465 




0.050 


0028125 


0010664 


0010159 






6p 


0.025 


0.0042014 


0.0041493 


0.0041548 






6d 


0.025 


0.0042014 


0.0040452 


0.0040606 






6f 


0.025 


0.0042014 


0.0038901 


0.0039168 






6g 


0.025 


0.0042014 


0.0036943 


0.0037201 







TABLE II: Energy eigenvalues of the Hulthen potential as a function of the screening parameter 
for 4p, 4d, 4f, 5p, 5d, 5f ,5g, 6p, 6d, 6f and 6g states in atomic units (H = m = e = 1) and for 
Z = l. 
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0.2 0.4 0.6 0.8 1 1.2 1.4 1.6 

5r 



FIG. 1: The variation of the effective Hulthen V e ff(r) and the approximate effective Hulthen 
potentials V e ff(r) with respect to the 5r for various values of the screening parameter. The 
parameters are in atomic units (H = e = m = 1) and 5 change from 0.4 to 5.2 in steps of 0.4. 
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